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SUMMARY 

The presence of radomes and instruments that are sensitive to water 
films or ice formations in the nose section of all-weather aircraft and 
missiles necessitates a knowledge of the droplet impingement character- 
istics of bodies of revolution. Because it is possible to approximate 
many of these bodies with an ellipsoid of revolution, droplet trajec- 
tories about an ellipsoid of revolution with a fineness ratio of 5 were 
computed for incompressible axisymmetric air flow. From the computed 
droplet trajectories, the following impingement characteristics of the 
ellipsoid surface were obtained and are presented in terms of dimension- 
less parameters: (l) total rate of water impingement, (2) extent of 

droplet impingement zone, (3) distribution of impinging water, and (4) 
local rate of water impingement. 


INTRODUCTION 

All-weather aircraft and missiles frequently have instruments 
located in the nose section of the fuselage that are sensitive to im- 
pinging atmospheric water droplets and ice accretion. For example, it 
has been found that the operation of an aircraft radar system located 
in a nose or wing radome is affected by a layer of ice or water dis- 
tributed over the radome surface. Therefore, it is necessary to evalu- 
ate, for given flight conditions, the expected distribution of various 
sizes of impinging water droplets over the nose section of the aircraft 
or missile. In addition, problems such as those encountered in the per- 
formance of external armament during flight in icing conditions require 
the evaluation of droplet impingement on bodies of revolution in order 
to determine where ice will form. 

Although a large variety of body shapes are used for radomes, rocket 
pods, and bombs, the impingement calculations may be made for a body 
selected to approximate a large group of these practical shapes. A pro- 
late ellipsoid of revolution is a good approximation for many of these 
bodies, and it has the additional advantage of a flow field that is 
known exactly for incompressible, nonviscous flow. 
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The trajectories of atmospheric water drople'ts about a prolate 
ellipsoid of revolution with a fineness ratio of 5 moving at subsonic 
velocities were calculated with the aid of a differential analyzer at 
the MCA Lewis laboratory. From the computed trajectories, the rate, 
distribution, and surface extent of impinging water were obtained and 
are sxnnmarized in terms of dimensionless parameters in this report. 


SYMBOLS 


The 

A 

a 

B 

CD 

d 

*^med 


following symbols are used in this report: 
semimajor axis 
droplet radius, ft 
semiminor axis 

drag coefficient for droplets, dimensionless 
droplet diameter, microns 
volume-median droplet diameter, microns 
collection efficiency, dimensionless 


K 


^ed 

L 

Re 


Reo 


2p\^a^U 

inertia parameter, — , dimensionless 

inertia parameter based on vo lime -median droplet diameter, 
dimensionless 

major axis of ellipse, ft 

local Reynolds number with respect to droplet, 2apg^v/p, 
dimensionless 

free-stream Reynolds number with respect to droplet, 2apg^u/p, 
dimensionless 


^®0,med 


free-stream Reynolds number based on volume -median droplet 
diameter, dimensionless 


r,z cylindrical coordinates, ratio to major axis, dimensionless 

tq starting ordinate at z = -<x of droplet trajectory, ratio to 

major axis, dimensionless 

ro tan starting ordinate at z = -os of droplet trajectory tangent to 
' ellipsoid surface, ratio to major axis, dimensionless 
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S distance along surface of ellipsoid from forward stagnation 

point to point of droplet impingement, ratio to major axis, 
dimensionless 

Sm limit of impingement zone, ratio to major axis, dimensionless 

t time, sec 

U free-stream velocity, ft/sec; or flight speed, mph, when indi- 

cated 

u local air velocity, ratio to free-stream velocity 

V local droplet velocity, ratio to free-stream velocity 

V magnitude of local vector difference between velocity of drop- 

let and' velocity of air, ft/sec 

W rate of impingement of water, Ib/hr 

W„ total rate of impingement of water on surface of ellipsoid, 

Ib/hr 

Wp local rate of impingement of water, lb/(hr)(sq ft) 

w liquid-water content in cloud, g/cu m 

a 1/4 focal distance of ellipsoid 

3 local impingement efficiency, dimensionless 

€ eccentricity of ellipse defined by a meridian section of 

ellipsoid of revolution 

A, |i prolate-elliptic coordinates 

viscosity of air, slugs/ (ft) (sec) 

density of air, slugs/ cu ft 

density of water, slugs/cu ft 

T time scale, tU/L, dimensionless 

Subscripts: 

r radial component 

z axial component 
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ANALYSIS 

The flow field around an ellipsoid of revolution in a stream moving 
in the direction of its major axis is axisymmetric. Therefore, the 
droplet impingement on the elliptical sections of all meridian planes 
is the same. Thus, the droplet impingement distribution on an ellipsoid 
of revolution can be obtained from trajectories in the z,r plane ^fig. 
l). The coordinates z and r are dimensionless and expressed as ratio 
of actual distance to major axis length L. The dimensionless equations 
of motion of the droplet trajectories are of the same form as those 
derived in reference 1 and can be written 


dT 


CpRe I 

— K K - ’'z) 




dT 


1 , > 
K - Vr) 


where 


(la) 


(lb) 


^ 2 Pw^^ 

9 liL“ 


( 2 ) 


and T = tU/L, and all velocities are in the form of ratio of local 
velocity to free-stream velocity U. 


The Reynolds number Re can be obtained conveniently in terms of 
the free-stream Reynolds number 


Reo = 2apgU/ji 


(3) 


from the relation 



(4) 


The coefficient of drag Cp is a function of Reynolds number. The 
values of Cp corresponding to various values of Reynolds number are 
obtained from experimental drag data (ref. 2). 


The air velocity components for incompressible nonviscous flow 
about a prolate ellipsoid of revolution are obtained from the exact 
solution of Laplace's equation in prolate-elliptic coordinates (fig. 2) 
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given by Lamb (ref. 3). The details for obtaining the velocity com- 
ponents in the z,r plane from Lamb's potential function in prolate- 
elliptic coordinates are given in appendix A. The z and r com- 
ponents of the air velocity field can be expressed in the form 

uz = - C I In 

and 

Ur = - f Cr ( 


where 

Jm 




C s 


■The constant € is the eccentricity of the ellipse defined by the 
meridian section of the ellipsoid of revolution. For an ellipse with 

fineness ratio of 5, € = <^.96. Equations (5) and (6) were solved for 
several hundred points in the flow field with the use of electronic cal- 
culating machines employing punched cards. The values of the air veloc- 
ity components U 2 and Ur as functions of r and z are given in 
figure 3. Figure 3(a) gives Uz as a function of z for constant 
values of r, while figure 3(b) gives Ur as a function of r for 
constant z . 

Assumptions that are necessary to the solution of the problem are: 

( 1 ) At a large distance ahead of the ellipsoid (free-stream condi- 
tions), the droplets do not move with respect to the air. 
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(2) No gravitational force acts on the droplets. 

(3) The droplets are always spherical and do not change in size. 

The first two assumptions are valid for droplets smaller than 
drizzle or rain drops. The assumptions are usually also valid for fall 
ing rain drops, because the airplane velocity is usually much greater 
than the drop velocity caused, by gravitational force. Preliminary cal- 
culations have shown that the third assumption is valid for the order 
of accuracy usually required in the design of equipment for the pro- 
tection of aircraft. 


METHOD OF SOLUTION 

The differential equations of motion (l) are difficult to solve, 
because the values of the velocity components and the term containing 
the coefficient of drag depend on the position of the drop at each in- 
stant and, therefore, are not known until the trajectory is traced. The 
values of these quantities must be fed into the equation as the tra- 
jectory of a droplet is developed. This was accomplished by using a 
mechanical differential analyzer constructed at the NACA Lewis labora- 
tory for this purpose -(ref. l). The answers were obtained in the form 
of plots of droplet trajectories with respect to the ellipsoid. A 
typical group of droplet trajectories is shown in figure 4. From the 
droplet -trajectory plots were obtained the impingement characteristics 
of interest discussed in subsequent sections. 

The equations of motion (l) were solved for various values of the 
parameter 1/K between 0.1 and 90. For each value of the parameter 
1/K, a series of trajectories was computed for each of several values 
of free-stream Reynolds number Roq: 0, 128, 512, 1024, 4096, and 8192. 

In order that these dimensionless parameters have more physical signifi- 
cance in the following discussions, some typical combinations of K and 
Reo are presented in table I in terms of - the length and the velocity 
of the ellipsoid, the droplet size, and the flight pressure altitude and 
temperature . 

Before the integration of the equations of motion to obtain the 
trajectories could be performed with the differential analyzer, the 
initial velocity of the droplets had to be determined at the point se- 
lected as the starting position. In addition, since the starting posi- 
tion must be selected at a finite distance ahead of the ellipsoid, it 
was necessary to make a correction to this starting ordinate in order 
to obtain the corresponding starting ordinate rQ at z = -00. Pre- 
liminary tests showed that, from z = -00 to z = -2, the trajectories 
and streamlines were essentially coincident. Furthermore, between 
z = -2 and z = -1, the magnitude of the change in the air velocity 
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components with change in r was negligibly small (fig. 3) compared 
with unity on the air velocity scale used in solving the equations. 
Consequently, the changes in the droplet velocity and the ordinate r 
along a given trajectory between z = -2 and z = -1 were essentially 
Independent of the starting r value or of the value of Roq* There- 
fore, a trajectory was run between z = -2 and z = -1 for each value 
of i/K in order to obtain the correct starting conditions at z = -1. 
By this method it was possible to solve the trajectories starting at 
z = -1 so that, within the accuracy of the machine, this would be 
equivalent to starting the trajectory at an infinite distance ahead of 
the ellipsoid. 

In order to minimize errors due to changes in the length of the 
paper used on the air velocity input drums caused by humidity and tem- 
perature variations, the component velocities of the flow field shown 
in figure 3 were plotted on 20- by 30-inch sheets of glass tracing 
cloth. A fine grid was laid out on the tracing cloth with a plotting 
machine constructed for this purpose. The droplet trajectories (fig. 4) 
were plotted by the differential analyzer on sheets of acetate in order 
to minimize scale changes and damage due to handling. The r-ordinate 
of the trajectory plots was scaled to 4 times the z-ordinate in an 
effort to Improve the accviracy of determining the point of droplet im- 
pingement on the ellipsoid surface. With this distorted scale, the 
trajectories were plotted with respect to an ellipsoid section with a 
major axis of 30 inches and a minor axis of 24 inches. The accuracy 
with which the various droplet impingement characteristics could be 
obtained is discussed in the following section. 


RESULTS AND DISCUSSION 

A series of droplet trajectories about the, ellipsoid of revolution 
with a fineness ratio of 5 at zero angle of attack was computed for the 
various combinations of the dimensionless parameters K and Rsq. 

These data are summarized in figure 5, where the starting ordinate rQ 

of each trajectory is given as a function of the point of inipingement 
on the surface S. (S is the distance measiired along the surface from 
the forward stagnation point to the point of impingement; the relation 
between S and z and r is given in appendix B.) The dashed lines 
in figure 5 are the loci of the termini of the constant K curves. 

These loci were found to be the same, within the order of accuracy of 
the computations, for all values of Rcq, as can be seen by comparing 
figures 5(a) to (f). The dot-dashed curves among the constant K curves 
for Reo = 4096 and Rcq = 8192 (figs. 5(e) and (f)) were obtained by 
interpolation. From the data presented in this figure, the rate, the 
area, and the distribution of water-droplet impingement on the surface 
of the ellipsoid can be determined for given values of Reo and K. 
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Total Rate of Impingement of Water 

In flight through clouds composed of droplets of uniform size, the 
total amount of water in droplet form impinging on the ellipsoid is 
determined hy the amount of water contained in the volume within the 
surface formed by the tangent trajectories (fig. l). Therefore, the 
total rate of impingement of water (ib/hr) can be determined from the 
relation 


W„ = 0.33,rg^t^^wL2u (7) 

where 0,33 is a conversion factor, the flight speed U is in miles per 
hour, the liquid-water content w is in grams per cubic meter, and L 
is in feet. When constants are combined. 


2 

In this equation, rQ is a measure of the efficiency of catch, be- 

cause it is proportional to the collection efficiency (Ejjj = lOOrQ^^g^^), 

defined as the ratio of the actual amount of water intercepted by the 
ellipsoid to the total amount of water in droplet form contained in the 
voliome swept out by the ellipsoid. 


The value of tq tan ^ given combination of Rcq and K can 

be obtained from figure 5 by determining the value of rQ which corre- 
sponds to the maximum S for the constant K curve of interest. The 
values of rQ tan fa-H the dashed termini curves of figure 5. In 

order to facilitate interpolation and extrapolation, the data are re- 


plotted in the form of 
in figure 6, 


^ , as a function of 
0,tan 

Examination of figure 6 shows that r, 


K for constant 
2 


Re, 


0 


0,tan 


increases with 


increasing K but decreases with increasing ReQ. 


The accuracy of the determination of Tq tan very much dependent 

on the shape of the tangent trajectory in the vicinity of the surface of 

the ellipsoid. For example, for the combination of i/K = 1 and 

Reo = 4096 shown in figiore 4, the shape of the tangent trajectory and its 

neighboring trajectories is such that a small increase of the order of 
0.0007 in rQ above the true value of will result in a trajec- 

tory that definitely misses the ellipsoid. A similar slight decrease in 
the value of rQ will result in a trajectory that unmisteikably impinges 

on the surface. The tangent trajectory is, therefore, relatively easy 



NACA TN 3099 


9 


to determine for this case} and ro^tan can, therefore, he determined 
to an accviracy of the order of ±0.0003. However, for a combination of 
i/k and Reo that results in small valvies of ro,tan> such as shown 
in figure 4 for i/K = 30 and Reo = 512, the crowding together of the 
trajectories near the ellipsoid and the tendency of the trajectories to 
have the same shape as the ellipsoid surface result in a possible ques- 
tion as to the true value of rQ^tan* Trajectories with rQ as much 

as 30 percent smaller than the true value of rQ^tan appear to be 

tangent. Therefore, in order to avoid selecting a trajectory that is 
not the tangent when determining ro,tan> value of ro was in- 
creased by small increments until a trajectory definitely missed the 
ellipsoid. Thus, an upper limit to the re^gion of possible tangency was 
established and used as a guide when selecting the tangent trajectory. 
With this method for determining the tangent trajectory, the accuracy 
of ro tan in this Reo and K region is within ±0.0007 for values 
of rQ >0.01. For reported values of ro,tan < 0.01, the accuracy 
in determining the tangent trajectory is somewhat indefinite, but appears 
to be within ±0.001. 

The effect of body size on the value of for selected cloud 

droplet size and flight conditions is illustrated in figure 7. The cal- 
culated values given in figure 7 for ellipsoids with a fineness ratio 
of 5 and major axis lengths between 3 and 300 feet are for flight at 50, 
100, 300, or 500 miles per hour through uniform clouds composed of drop- 
lets of 10, 20, or 50 microns in diameter at pressure ’ altitudes of 5000, 
15,000, or 25,000 feet and temperatures (most probable icing temperature 
given in ref. 4) of 20°, 1°, and -25° F, respectively. For example, 
consider a 40-foot-long ellipsoid with a fineness ratio of 5 traveling 
at 500 miles per hour with zero angle of attack at a pressirre altitude 
of 15,000 feet through a uniform cloud composed of droplets, of 20 microns 
in diameter. From figure 7(b), Tq is 0.000114. If the liquid- 

water content of the cloud is assumed to be 0.1 gram per cubic meter, 
then (from eq. (S)) the total rate of impingement of water Wj^ is 9.5 

pdunds per hour. 


Extent of Droplet Impingement Zone 

The extent of the droplet impingement zone on the surface of the 
ellipsoid is obtained from the tangent trajectories. The point of tan- 
gency detennines the rearward limit of the impingement zone. The limit 
of impingement for a particular Roq and K condition can be 

determined from the maximum S value of the constant K curve of 
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interest in figure 5. Again, to facilitate interpolation, the data are 
replotted in the form of Sm as a function of K for constant Rcq 

values in figure 8. The data of this figure indicate that in- 

creases with increasing K but decreases with increasing ReQ- 

Because of the difficulty of determining the exact point of tan- 
gency on the surface of the ellipsoid of each tangent trajectory, the 
accuracy of determining Sjjj is of the order of ±0.005. The accuracy 
of determining the value of S for the intermediate points of impinge- 
ment given in figiire 5 was much higher, because the points at which the 
intermediate trajectories terminated on the ellipsoid surface were much 
better defined. 

The effect of body size on the value of for selected cloud 

droplet and flight conditions is illustrated in figure 9. For example, 
consider a 20-foot-long ellipsoid with a fineness ratio of 5 traveling 
300. miles per hour at zero angle of attack at a 5000-foot pressure 
altitude through a uniform cloud composed of 50-micron droplets. From 
figure 9(a), is 0.109; that is, the impingement zone extends 2.18 
feet rearward (measured along the surface) from the forward stagnation 
point. 


Distribution of Impinging Water Along Ellipsoid Surface 

The amount of water impinging on the ellipsoid surface within any 
ring of width S 2 - S]_ can be determined if the starting ordinates rQ 
are known for the droplets that impinge at Sj_ and Sg. These data 
can be obtained from figure 5. 

The amount of water (ib/hr) impinging within the ring of width 
S 2 - is given by the relation 

W = 1.04(r2^2 - O) 


Local Rate of Impingement of Water 


The local rate of impingement of water in droplet form 
(lb/(hr)(sq ft)) on the surface of the ellipsoid can be determined from 
the expression 


>■0 '^0 

Wp = 0.33U» - gg- 


= 0.33Uw(3 


(lOl 
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where p is the local impingement efficiency. The values of P as a 
function of S for combinations of Reo and K are presented in fig- 
ure 10. These c\irves were obtained by multiplying the slope of the 
curves in figure 5 by the corresponding ratio tq/t at each point. The 
dot-dashed curves included for Roq = 4096 and Roq = 8192 were obtained 
from the corresponding interpolated curves of figure 5. Because the 
slopes of the tq against S curves (fig. 5) in the region between 
S = 0 and S = 0.01 are difficult to determine, the exact values of P 
between S = 0 and S = 0.01 are not known. The values of P presented 
in this region between S = 0 and S = 0.01 are estimated to be accurate 
within iO.05. 


IMPINGEMENT IN CLOUDS OF NONUNIFORM DROPLET SIZE 


The data presented in figures 5 to 10 would apply directly only to 
flights in clouds composed of droplets that are all uniform in size. 

The droplets in a cloud, however, may have a range of sizes. Theoreti- 
cal calculations (ref. 5) and experience in the NACA Lewis icing re- 
search tunnel on bodies of revolution have shown that the amount of ice 
collected when a distribution of droplet sizes is present in the tunnel 
is considerably greater than that which would be obtained if only drop- 
lets of the volume -median size were present. Therefore, if the cloud 
droplet-size distribution is known or can be estimated, then the data 
must be accordingly modified (or weighted) before the rate, the extent, 
and the distribution of droplet impingement on the ellipsoid are 
calculated. 


For a nonuniform, cloud, the total rate of Impingement of water on 
the ellipsoid can be determined from equation (8) by using the weighted 
value of r2 that corresponds to the droplet-size distribution 


‘’0,tan 

present in the cloud. The weighted value of r,^ 
by plotting 


^ . can be obtained 
0,tan 

Tq for each droplet size (based on values of 


K and 

ReQ corresponding to each droplet dlam. ) as a function of the cumula- 
tive volume (in percent) of water corresponding to each droplet size 
and Integrating the resultant curve. For example, consider the cloud 
droplet-size distribution shown in figure 11. Suppose that the volume- 
median droplet size is 20 microns, the velocity is 200 miles per hour, 
the ellipsoid length is 10 feet, the pressvire altitude is 5000 feet, and 
the temperature is 20° F. For these conditions, the value of ReQ^med 

is 117.6 and of is 0.03898. The values of Roq and K corre- 

sponding to other droplet sizes in the distribution are obtained by 
multiplying Reo^med ^ed (V<^ed)^ 
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to ottain rg (fig. 6) for each droplet size. The values of 

7-2 for this example are plotted as a function of cumulative volume 

0,tan 

'(in percent) in figure 12. Integration of this curve gives a weighted 

_ O r\r\r\Anzz^ t rl-. c* +Vici trial no nn i".VlP 


value of 


O.tan 


equal to 0.000425; whereas, the value "based on the 


volume-median droplet size is 0.00031 (fig. 6) 


local rate of impingement of water at any point for a distrihu- 
tion of droplet sizes can he obtained in the same manner. The extent 
of the droplet impingement zone should he determined from the values of 
iC and ReQ calctaated for the largest droplets present in sufficient 

number to represent a significant portion of the total water present 
in the cloud. 


CONCUJDING REMARKS 

The scale factors used in the differential analyzer to solve the 
equations for the range of conditions presented in figures 5 to 10 and 
the near-parallelism of the trajectories to the surface at large values 
of 1/K made it impossible to obtain sufficient accuracy to present 
detailed data, such as the rate of local Impingement of water, at points 
along the surface of the ellipsoid, for values of 1/K > 90 for Roq = 0 

and 1/K > 30 for Re^ > 128. From table I it can be seen that, for 

bodies as large as the fuselage of cargo or passenger airplanes, these 
conditions are not uncommon. Examination of figures 5 and 10 shows, 
however, that the extent (usually < 0.03) and rate of local impinge- 
ment are small in this ReQ and K region. Therefore, in this region 
a knowledge of the extent of the impingement zone and the total rate of 
ing)ingement of water as calculated from the data of figures 8 and 6, 
respectively, is sufficient for most applications. 

Because the droplet trajectories about the ellipsoid were calcu- 
lated for incompressible fluid flow, a question as to their applicability 
at the higher subsonic velocities may arise. In reference 6 it was 
shown that the effect of compressibility up to the flight critical Mach 
number on the trajectories about a cylinder was negligible. In view of 
the results obtained for the cylinder and of the high flight critical 
Mach number (greater than 0.9) for the ellipsoid, the ellipsoid impinge- 
ment results should be applicable for most engineering uses throughout 
the subsonic region. 

The data of this report apply directly only to ellipsoids of revo- 
lution with a fineness ratio of 5. Therefore, consideration must be 
given to the degree of geometric and aerodynamic similarity before 
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applying the data to "bodies of revolution with other shapes. In some 
cases, where the body is of different shape, it may be possible to match 
its nose section physically with the nose section of an ellipsoid (fine- 
ness ratio, 5) of selected length. If, in such a case, the contribution 
of the afterbody to the air-flow field in the vicinity of the nose of 
the body is small (as. it often is), then the impingement data for the 
matching portion of the surface of the ellipsoid can be used for deter- 
mining the impingement characteristics of the nose region of the body. 

In other cases, where the body shape differs from that of an ellipsoid 
but the fineness ratio is the same, the air-flow field may be similar 
enough that an estimate of the total catch can be obtained from the 
ellipsoid data. In this case, no details of the surface distribution 
of impinging water could be obtained. . 


Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, November 20, 1953 
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APPENDIX A 

CALCULATION OF VELOCITY FIELD ABOUT AN ELLIPSOID OF BEVOLUTION 

The incompressible nonviscous velocity field for axisymmetric flow 
about an ellipsoid of revolution can be obtained as follows: 

Consider a prolate-elliptic coordinate system in the z,r plane 
(fig. l) defined by 


2a , t 

z = — cosh ^ cos TJ 

1j 


2a . , „ 

r = -=- sxnh C sin t] 
Li 




(Al) 


where 


0 < C 5 * 0 < T) < 2it 

The coordinates z and r are dimensionless and are expressed as 
a ratio to the major axis of the ellipsoid of revolution of interest. 
Examination of equation (Al) shows that ^ = constant and t| = constant 
represent confocal ellipses and hyperbolas, respectively, with a dimen- 
sionless semifocal distance of 2a/L. 


Let 


X = cosh ^ 


and 


H = cos T] 


Then 





r 




and 


> 


(A2) 


NACA TN 3099 


15 


where 


1 < X 00 and -1 ^ 1^ +1 

and X = constant and \x = constant also represent confocal ellipses 
and hyperbolas, respectively (fig. 2). The major axis of each ellipse 
is along the z-axis. 

The velocity potential expressed in X,|i coordinates can he ob- 
tained from zonal siirface harmonics (refs. 3, 7, and 8) and expressed 
in the following dimensionless form: 


where Viflg the usual dimensions of velocity potential. 

For an ellipsoid defined by the surface Xq = constant in a fluid 

moving with a free-stream velocity of U in the direction of the z-axis 
of the ellipsoid, the constant C is given by 


2a 


^0 1 , 

X2-l " 2 >.0-1 


(A4) 


The coordinates X and p can be e:q)ressed in terms of z,r coordinates 
as follows : 




(A5) 


and 




(A6) 


where C = ^ is the eccentricity of the ellipse defined by 

Xq = constant. If the surface coordinates of the ellipse of eccentricity 

€ are substituted in (A5), it is found that 


^0 ~ € 
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Then equation (A3) can he expressed in the form 


^ (I ^ In ^ - l) 



where 





(A7) 


In addition, equations (A2) take the form 


z = 5. 
2 


and 


= |a/x2 - 1^1 - h2 


The A and p. velocity components are obtained from the relations 


^A = 




2 A x2-i 


5a 


2 /A^-l 


'A2.^x2 


1^2 x:t 


A^-l 


2 


and 


(A8) 






-1 




£ IS! 
2 


2 /]^ 
€ I 


a2-^.2 


^ , 1 -N ^ 7\+l 

C I 2 Mti xir 




(A9) 
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The z and r velocity components in the z,r coordinate system 
are obtained from equations (A8) and (A9) as follows: 



(All) 


The velocity components and Uj. can he written as a function 

of z and r hy substituting equations (A5) and (A6) into equations 
(AlO) and (All): 



and 


Ur 



1 .] i 


( 6 ) 


where 




and C is given in equation (A7). 
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The solution of equations (s) and (6) at several hundred points in 
the flow field for an ellipsoid of revolution with a fineness ratio of 

5 (e = //07^) was accomplished with the use of electronic calculating 
machines employing punched cards. The values of u^ and Up as func- 
tions of r and z are given in figure 3. Figure 3(a) gives Uz as 
a function of z for constant values of r, and figure 3(h) gives uj. 
as a function of r for constant z. 
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APPENDIX B 


RELATION BETWEEN DISTANCE ALONG SURFACE OF ELLIPSE AND z AND r 

The length of arc of an ellipse cannot be reduced to an elementary- 
function of z and r, but rather belongs to the class of functions 
known as elliptic integrals. The length of arc S as a function of 
z and r is obtained (by the method of ref. 9) as follows: 

The equation for an ellipse in the z^r plane can be written 
parametrically in terms of d as follows: 


z = A sin 6 I 

f (Bl) 

r = B cos 0 1 

where A and B are^ respectively, the semimajor and semiminor axes 
of the ellipse (fig. l). 


Then 


dS^ = dz^ + dr^ 



cos^0 + B^ 


sin% 


d0' 



(A^-B^) sin^e 


d02 


and 


S 


e 



(B2) 


where e 


Va2-b2 

A 


denotes the eccentricity. 


This function is of the form 




1 - sin^G d0 


(B3) 


where 


0 < k < 1 


which is known as the Elliptic Integral of the Second Kind in Legendre ' s 
form. Tabulated values of this function are available (ref. lO). 
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For the purposes of this report, the distance along the surface 
measured from the nose of the body (z = -0.5 and r = O) is desired, 
or. 


S = A - k2 sin^e d0 




(B4) 


where 


A = 0.5 


and 


^2 _ g2 _ 0.96 

The relation between S and z computed from equation (B4) is 
shown in figure 13. The relation between S and r^ is given by the 
1/K = 0 curve of figure 5, inasmuch as r = rQ for infinitely large 

values of K. 
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Figure 2. - Prolate-elliptic coordinate system. 
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(a) z-Component of air velocity as function of z for constant values of 
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(b) r-Component of air velocity as function of r for constant values of 
Figure 3. - Concluded. Ellipsoid flow field. 
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Figure 4. - Trajectories of droplets with respect to ellipsoid. 
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Figure 5. - Starting ordinate as function of distance along surface to point of Impingement. 
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(b) Pree-stream Reynolds number, 128. 

Figure 5. - Continued. Starting ordinate as function of distance along surface to point of Impingement. 
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Figure 5. - Continued. Starting ordinate as function of distance along surface to point of Impingement. 
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Figure 5. - Continued. Starting ordinate as function of distance along surface to point of Impingement 
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Figure 6. - Square of starting ordinate of tangent trajectory as function of Inertia parameter. 
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Figure 7. - Square of starting ordinate of tangent trajectory as function of major 
axis length of ellipsoid. 
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(b) Pressure altitude, 15,CX)0 feet; temperature, 1° F. 


Figure 7. - Continued. Square of starting ordinate of tangent trajectory as function 
of major axis length of ellipsoid. 
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(c) PreaBure altitude, 25,000 feet; temperature, -25° P. 


Figure 7. - Concluded. Square of starting ordinate of tangent trajectory as function 
of major axis length of ellipsoid. 
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Figure 8, - Limit of impingement zone on surface of ellipsoid 
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Figure 9. - Maximum extent of Impingement zone as function of major axis of ellipsoid. 
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Figure 9. - Continued. Maximum extent of Impingement zone as function of major axis 
of ellipsoid. 
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Figure 9. - Concluded. Maximum extent of Impingement zone as function of major axis 
of ellipsoid. 
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Figure lO. - Local Impingement efficiency; 
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Figure 10. - Continued. Local Impingement efficiency; 
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Figure lO. - Continued. Local Impingement efficiency: 


MCA TN 3099 


45 



O <7> 


O 


d *^ouaf ofjja i^uauiaSufduJT Iboot 


.16 .20 .24 .28 .32 .36 .40 .44 .48 .52 .56 

Distance along surface to point of Impingement, 3 


46 


NACA TN 3099 




II 

(H 


Figure 10. - Continued. Local Impingement efficiency; 
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Figure 10. - Concluded. Local impingement efficiency: 
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Figure 11. - Typical cloud droplet-size distribution In NACA Lewis Icing 
research tiumel. Vol\ime-medlan droplet diameter, 20 microns. 
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Figure 12. - Relation of to cumulative volume of water corresponding to droplet size from 

which To,tan determined. Area under curve gives (r^ " 0.000425. 
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